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A novel and robust pedestrian detection method in thermal infrared images based on the double-density 
dual-tree complex wavelet transform (DD-DT CWT) and wavelet entropy is presented in this paper. The 
regions of interest (ROIs) are located first making use of high brightness property of the pedestrian pixels 
caused by the self-emission of the pedestrians related to the Planck’s law. The candidate ROIs are then 
decomposed by DD-DT CWT and the wavelet entropy features are extracted from the high frequency sub¬ 
bands. The true pedestrian regions are finally classified and recognized using the support vector machine 
(SVM) classifier. Comparisons between our approach and traditional approaches are presented and 
experimental results using several thermal infrared image databases show the proposed scheme to be 
very promising. 

© 2010 Elsevier B.V. All rights reserved. 


1. Introduction 

Recently, automatic pedestrian detection in thermal infrared 
imagery for intelligent video surveillance and driver assistance- 
systems for vehicles has attracted more and more attention [1- 
4]. However, the problem of detecting pedestrians is very difficult 
because of the variability in appearance and pose for each pedes¬ 
trian, especially in outdoor environments. Compared with visible 
light images, thermal infrared images have unique characteristics. 
Generally, the target object’s intensity is mainly determined by its 
temperature and radiated heat and is independent of the current 
light conditions, so the detection system can be applied indiscrim¬ 
inately in both day and night. Also, the infrared images almost 
eliminate the influences of color, texture, and illumination on the 
target object’s appearance variability. Therefore, there is a latent 
development potential for pedestrian detection in thermal infrared 
images. However, infrared images are not perfect either. First, be¬ 
cause the target object’s surface properties (emissivity, reflectivity, 
and transmissivity) and wavelength both affect the infrared images 
intensities, non-human objects and backgrounds, such as animals, 
cars, light poles, and buildings, produce additional bright areas in 
infrared images, especially in summer afternoon. These clutters 
in infrared images make it impossible to accurately detect pedes¬ 
trians based only on their brightness. Second, due to limitations 
in camera technology, most infrared images have lower spatial res¬ 
olution and less sensitivity than visible images which often lead to 
low image quality, such as blurring, low target-to-background con¬ 
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trast, and great noise. So it is a complex challenge to make precise 
detection of pedestrians using thermal infrared imagery. 

In general, pedestrian detection in thermal infrared images is 
performed in two steps: candidate ROI extraction and pedestrian 
detection from these candidate regions. The main methods in the 
first step include template match [5], movement difference [2,6] 
and intensity oriented projection [7-9]. The template match meth¬ 
od requires multiscale searching for the whole image, which is 
time-consuming. The movement difference method is mainly ap¬ 
plied to moving images and is not appropriate for static or single 
images. This paper therefore uses a statistical adaptive intensity 
oriented projection method to extract regions of interest. The sec¬ 
ond step adopts certain pattern recognition or machine learning 
methods to classify the candidate ROIs and recognize the real 
pedestrian regions. The key problem is how to describe the pedes¬ 
trian and discriminate the pedestrians from other candidates based 
on effective features. Many features are currently used to describe 
pedestrians, such as shape features (compactness and leanness) 
[10], histograms of oriented gradients (HOG) [11,12] and shape- 
independent features (histogram, inertia, and contrast) [9]. This 
paper uses wavelet entropy features based on the double-density 
dual-tree complex wavelet transform (DD-DT CWT) to accurately 
describe the pedestrian. 

In this paper, we introduce a complete robust pedestrian detec¬ 
tion system for application in thermal infrared images. Based on 
the self-emission of pedestrians related to the Planck’s law, which 
usually results in the high brightness property of the pedestrian 
pixels in thermal infrared images, we first adopt the intensity ori¬ 
ented projection method to extract the candidate ROIs which 
avoids the time-consuming multiscale searching of the whole 
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image. We then acquire the wavelet entropy features of candidate 
ROIs by using DD-DT CWT, which can precisely describe the pedes¬ 
trian’s features. Finally, we use the SVM classifier to recognize the 
true pedestrians from the candidate regions. Detection results are 
found for the OTCBVS Benchmark Dataset Collection [13] and com¬ 
parisons between our approach and traditional detection methods 
are presented. Experimental results show that the proposed 
scheme can improve the robustness of the complete system and 
precisely detect the pedestrians. 

2. Candidate pedestrian regions selection 

During the selection of candidate pedestrian regions, we use a 
statistical adaptive intensity oriented projection method which 
contains horizontal and vertical projections. In order to eliminate 
the interference from noise and low brightness pixels, we must 
first segment the original image into a binary image using a flexible 
threshold T. Considering that the dynamic range of image intensity 
between frames may vary dramatically due to different seasons or 
varying outdoor environments, we set the threshold T as a balance 
between the image maximum intensity I max , the image mean 
intensity I mean and the pedestrian mean intensity P mea n in the train¬ 
ing sample set, expressed as 

T = Wi I max + W 2 Imean + (1 - W t - W 2 )P mean (1) 

where Wi and w 2 are weights that satisfy 0 < wi, w 2 < 1 . The values 
of Wi and w 2 are set to acquire the optimal segmentation results 
through experimental tests. The threshold T will vary adaptively 
with the dynamic range of image intensity. Then the segmented 
binary image is projected vertically using the following steps: 

(1) Obtain the intensity vertical projection curves by computing 
a count of the pixels whose intensities are higher than the 
threshold T. The intensity vertical projection curves are 
defined as the number of bright pixels in image columns 
according to their horizontal positions, as shown in Fig. lb. 
Fig. la is the original image. The curves can be divided into 
several heaves or waves with rising left curves and declining 
right curves. 

(2) Search for each rising left edge and declining right edge from 
the projection curves. 


(3) Segment the image into several vertical stripes by pairing 
each rising left edge and declining right edge, as shown in 
Fig. lc. Each stripe may contain one or more pedestrians. 

Like the intensity vertical projection curve, the horizontal pro¬ 
jection curve can also be defined as the number of bright pixels 
in image rows according to their vertical positions. A horizontal 
projection is used to acquire horizontal stripes based on the projec¬ 
tion curve. Then we can make the bounding boxes from the inter¬ 
section of the vertical and horizontal stripes. These bounding boxes 
are called the candidate ROIs where pedestrians may be located. 

However, the candidate ROIs may contain some unwanted arte¬ 
facts due to interference from the many warm targets in the scene. 
Therefore, an additional intensity oriented projection method is 
applied to produce more accurate candidate pedestrian ROIs using 
the following steps: 

(1) Compute the width/height ratio of the candidate ROI as 
Rwh = Uk, where l w and l h represent the width and height 
of the candidate ROI. 

(2) If R wh ^ T v , then the candidate ROI is projected by the sec¬ 
ond level of vertical oriented projection, where T v represents 
the second level of vertical projection threshold. 

(3) Else if R wh < T h , then the candidate ROI is projected by the 
second level of horizontal oriented projection, where T h rep¬ 
resents the second level of horizontal projection threshold. 

Fig. 1 shows the results of the intensity oriented projection 
method for a thermal infrared image using two levels of projec¬ 
tions. We can see from Fig. Id and e that the second level of projec¬ 
tion almost eliminates the interference from the trailer 
phenomenon and acquires more accurate candidate ROIs. 

3. Feature extraction 

Feature extraction is the key point of the complete pedestrian 
detection system. The wavelet transform is often used for robust 
feature extraction from images or time series data for its advantage 
of local analysis in the spatial and frequency domains. Recently, the 
dual-tree complex wavelet transform and double-density wavelet 
transform were both used to extract the mean and variance of 





Fig. 1 . Results of the intensity oriented projection method for a thermal infrared image: (a) original image, (b) intensity vertical projection curves, (c) horizontal 
segmentation, (d) first level of projection, and (e) second level of projection. 
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the wavelet coefficients from each scale and frequency band, and 
were successfully applied to palmprint recognition [14], texture 
classification [15], and image retrieval [16,17]. This paper proposes 
a novel feature extraction method based on DD-DT CWT, which 
combines both the dual-tree wavelet transform and the double¬ 
density wavelet transform. DD-DT CWT describes the pedestrian 
image with 16 dominant orientations and each orientation is rep¬ 
resented using two wavelets, therefore it can improve the accuracy 
of image decomposition and reconstruction. Also, the wavelet en¬ 
tropy is taken as the feature vector by computing the wavelet coef¬ 
ficients from each transformed frequency band, which reflects the 
energy distribution in the frequency domain. The final wavelet en¬ 
tropy feature thus reduces the dimensions greatly and therefore 
improves the speed of the detection system. 

3.1. Double-density dual-tree complex wavelet transform 

The 2-D dual-tree complex wavelet transform (2-D DT CWT) 
[18,19] is implemented by using four 2-D discrete wavelet trans¬ 
forms in parallel with different filter banks in rows and columns 
separately. Twelve wavelets are acquired by taking the sum and 
difference of each pair of subbands. Six directions are described, 
each by two wavelets. 

The 1-D double-density wavelet transform (1-D DD WT) [20] is 
based on a scaling function and two different wavelets where one 
wavelet is a half-sample shift of the other. Compared with tradi¬ 
tional wavelet transforms, oversampling is used instead of critical 
sampling during implementation to reduce the translation sensi¬ 
tivity. The 2-D DD WT [21] uses the 1-D DD WT alternately; the 
row is filtered first, and then the column. Nine 2-D subbands are 
created, where one is a 2-D subband and the remaining eight sub¬ 
bands become eight 2-D wavelet filters. The orientation descrip¬ 
tions are like those of the discrete wavelet transform. 

The 1-D DD-DT WT [22] combines the 1-D DT WT and the 1-D 
DD WT. It is constructed using two different scaling functions 
(</>h(t), 4> g ) and four different wavelets (i \j h ^ ij/gj, i =1,2), where 
is a half-sample shift of i/q, j2 (t) and ij/ g ^{t) is a half-sample 
shift of i/^ i2 (t), as follows: 

'I'hAh ='I'h , 2 (t — ° 5 ), ^,i(t)*^(t-0.5) (2) 

The two wavelets become an approximate Hilbert transform 
pair as 

= K 2(t) = H{iMt)} (3) 

The 2-D DD-DT CWT is implemented using four oversampling 
2-D DD DWTs for the input image in parallel. The rows and col¬ 
umns are filtered using different filter banks. Fig. 2 illustrates the 
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process of two levels of the 2-D DD-DT CWT. Lo p and Hi p are the 
first level of the decomposition filter banks and represent one scal¬ 
ing filter and eight wavelet filters of 2-D DD DWT. Lo and Hi are the 
second and remaining levels of decomposition filter banks. Lo mn 
and Hi mn (m = 1,2; n = 1,2,..., 4) represent one low pass subband 
and eight high pass subbands separately. Therefore, four low pass 
subbands and 32 high pass subbands are produced for our input 
image after one level of the transform. 

The DD-DT CWT is completed by iteratively transforming the 
low pass subbands with the determined levels. Finally, 32 wave¬ 
lets are created by the sum/difference operation on each pair of 
subbands. Sixteen orientations are described, each by two 
wavelets. 

Fig. 3 shows the comparisons between the 2-D DT CWT, 2-D DD 
WT and 2-D DD-DT CWT. In Fig. 3a, six directions (±15°, ±45°, ±75°) 
are described by 2-D DT CWT using twelve wavelets in the first and 
second line, each by two wavelets (one for real part, the other for 
imagery part), and the third line can be taken as the magnitude 
of the complex wavelet. Fig. 3b illustrates eight wavelets used in 
2-D DD WT, the first two wavelets are oriented in the vertical 
direction, the third and sixth are oriented in the horizontal direc¬ 
tion, while the remaining four wavelets have no dominating orien¬ 
tation for the checkerboard affect. Fig. 3c shows 32 wavelets used 
in 2-D DD-DT CWT, which describe 16 orientations. As we can see, 
after the 2-D DD-DT CWT, the frequency band in the transform do¬ 
main is finer, and the orientation information is richer. Also, the 2- 
D DD-DT CWT eliminates the mosaic phenomenon and acquires 
more detailed information for the edges of the input image. There¬ 
fore we can get more accurate representation for pedestrian re¬ 
gions in thermal infrared images after 2-D DD-DT CWT. 

The candidate pedestrian regions are not large in the thermal 
infrared images, so we set the wavelet decomposition scale factor 
j = 3. Suppose the candidate region size is mx n; we then get 
three levels of subbands whose sizes are m/2 J x n/2 J . Thirty-two 
high pass subbands and four low pass subbands are produced 
after each level of decomposition. Three levels of decompositions 
are completed by iteratively transforming the low pass subbands. 
Finally, 96 high pass subbands and four low pass subbands are 
created in total. Fig. 4 shows the decomposition processes of a pe¬ 
destrian and car in thermal infrared images using 2-D DD-DT 
CWT. Two subbands are used to describe the detail information 
for the same dominant orientation in each level and therefore 
there are 16 main orientations in all. Each high frequency sub¬ 
band represents the fine-to-coarse description from different ori¬ 
entations and resolutions on the same side of the edges and 
contours. 



Fig. 2. Two levels of 2D DD-DT CWT. 


Fig. 3. (a) 2-D DT CWT, (b) 2-D DD WT, and (c) 2-D DD-DT CWT. 





















































































270 


J. Li et al. I Infrared Physics & Technology 53 (2010) 267-273 


a 



b 



Tree One 


IajvxjI One /^jj| 

r r 

Level Two »j 

UvelTh'ee », *, 

Frequency * 


Tree Two 



i) i j i 


|/| 

[ t 


Tree One 


Level One 


Level Two 


Level Thax 
Low 

Frequency 




Tree Two 



I I 


Fig. 4. Decomposition processes of a pedestrian and car using 2-D DD-DT CWT: (a) pedestrian decomposition; (b) car decomposition. 


3.2. Wavelet entropy 


4. SVM classifier 


Wavelet entropy is widely used as a feature extraction method 
in signal detection and recognition applications. Let the wavelet 
energy Ej be the quadratic sum of the wavelet coefficients Cj{k) in 
a single scale j as 

Ej = jr\C i (k)\ 2 (4) 

k =l 

The total energy of the signal is 

£**/ = £E, = ££lQ(k)| 2 (5) 

j =1 J =1 k=i 


The support vector machine (SVM) classifier is a binary classi¬ 
fier algorithm that tries to find an optimal hyperplane as a decision 
function in a high dimensional space [24]. It has great generaliza¬ 
tion ability and is adept at solving pattern recognition problems for 
non-linear and small sample sets. In this paper, we use a SVM clas¬ 
sifier to classify and recognize pedestrians. 

In a two-class pattern recognition problem, given a training 
data set: 

{(XjJ i );X i eR N ;y i = ±l;i=l,2.J} (9) 

the decision function is set as 


The relative wavelet energy is defined as 


pj 



EPj = 1 


j= 1 


(6) 


If the wavelet coefficients matrix is taken as a probability distribu¬ 
tion series, then the entropy computed by the series reflects the 
sparseness of the coefficients matrix, which is the order degree of 
the distribution derived by wavelet entropy [23]. According to the 
Shannon entropy theory, we define the wavelet entropy as 

S w = -'£ p jlnP j (7) 

3=1 


The wavelet entropy of the pedestrian image reflects the energy dis¬ 
tribution of the image in the frequency domain of the wavelet 
transform. It is defined as 

J 

H m (j) = (8) 

3=1 

where j is the scale, and m represents the high frequency part. In the 
2-D DD-DT CWT, m = 1,2 ... 32, which includes 16 directions of the 
real part and 16 directions of the imaginary part. A 1 x 32 vector is 
thus produced by extracting the wavelet entropy from the original 
pedestrian image after 2-D DD-DT CWT, which greatly reduces the 
feature dimensions. 


/(x) = sgn 


ya,K(xi,x) + b 


( 10 ) 


_i= 1 

with the constraint conditions given as: 

! i 

min - (co, co) + C C ^ 0; 

co( 2 tf 

y,[(co,x) + b] > 1 {,• » 0. (11) 

The optimal solution of this problem is then obtained using the 
Lagrangian theory 

i j 

max + o jK(X[,Xj), 0 < a* < C (12) 

i-\ 1 


where x, is the training sample set, y t is the class label, x is the pend¬ 
ing sample, and b is the threshold. I<( •) is the kernel function which 
satisfies the Mercer condition. The role of I<( •) is to map the original 
feature space to high dimensional space and even to infinite dimen¬ 
sional space. The optimal solution of the problem is then sought in 
the new space. Popular kernel functions include the Gauss kernel, 
the polynomial kernel, and the Sigmoid kernel. 

Training the SVM classifier involves looking for the solutions for 
a and b that maximize the formula (12). This optimization problem 
can be tuned by changing the value of the constant C and the 
parameters of the kernel function. The Gauss kernel function is 
chosen and the parameter optimized method proposed by Carl 
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Fig. 5. Block diagram of detection system. 
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Table 1 

Experimental results of different features based pedestrian detection algorithms. 



HOG 


HOG + SHAPE 


HOG + INERTIA 

HOG + SHAPE + INERTIA 

WE (DD-DT CWT) 

FPR (%) 

TPR (%) 

FPR (%) 

TPR (%) 

FPR (%) 

TPR (%) 

FPR (%) 

TPR (%) 

FPR (%) TPR (%) 

1 

12.6 

91.0 

3.9 

94.1 

7.9 

94.1 

3.9 

98.1 

1.6 98.7 

2 

1.7 

94.4 

1.1 

95.8 

1.8 

94.1 

1.1 

96.4 

0 99.2 

3 

1.3 

94.1 

2.7 

95.1 

1.1 

97.8 

1.0 

98.0 

0.05 99.0 

All 

1.9 

93.9 

2.1 

95.4 

1.6 

95.8 

1.2 

97.2 

0.09 99.1 



















0995 





099 

. cDD-DT CWT 


















099 

(DO-OTCWT 













<DT CWT 



0.98 







0985 

dDOCrwr 







(WOG+SHAPE+INERTIA 




096 





a, 0.97 







fe 





•5 







% 





at 







a 0 975 





a. 







Q. 





Y- 0 96 







H 








<WOG+ INERTIA 



097 









(WOG+SHAPE 







095 







0965 















(DWT 









096 





094 




<WOG 








093 

_i_ 

_i_ 

_i_ 

_i_ 



0956 

°*< 





( 

) 0005 

0.01 

0.015 

0.02 

0.025 

1 0005 0 01 0.015 

0.02 0.025 

0.03 0035 0.04 0.045 005 


FP Rale 

Fig. 6. Comparisons between different algorithms. 


FP Rate 

Fig. 9. Comparisons between different wavelet transforms. 




and Peter [25] is used in this paper, and therefore each parameter 
is tuned adaptively in order to obtain the optimal classification 
results. 

5. Complete system 

In this section, we describe a proposal for a complete pedestrian 
detection system for thermal infrared images. The framework of 
the system is shown in Fig. 5, which contains three models: the 
candidate regions selection, feature extraction and classification. 
The entire procedure is listed in the following steps: 

(1) Search for all possible candidate ROIs from the input image 
using the two-level intensity oriented projection method. 


(2) Classify the candidate regions manually, and randomly 
select 325 pedestrian samples and 336 non-pedestrian sam¬ 
ples which are used to construct the training data set. 

(3) Decompose all the training samples by three levels of 2-D 
DD-DT CWT and therefore obtain 96 high frequency sub¬ 
bands and four low frequency subbands from each training 
sample. 

(4) Acquire the wavelet entropy by computing the relative 
energy of each high frequency subband, where a 1 x 32 vec¬ 
tor is created from each training sample. 

(5) Take the 1 x 32 vector acquired in step 4 as the input feature 
of the SVM classifier, adaptively tune all of the parameters of 
the classifier and save the SVM classifier with the best recog¬ 
nition results for the remaining detection process. 
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Table 2 

Experimental results of different wavelet transforms. 



DWT 


DT CWT 


DD DWT 


DD-DT CWT 


FPR (%) 

TPR (%) 

FPR (%) 

TPR (%) 

FPR (%) 

TPR (%) 

FPR (%) 

TPR (%) 

1 

7.9 

92.2 

3.1 

96.1 

3.9 

96.1 

1.6 

98.7 

2 

6.1 

99.0 

2.9 

99.2 

3.8 

99.0 

0 

99.2 

3 

3.7 

94.1 

1.6 

98.8 

0.8 

98.6 

0.05 

99.0 

All 

4.8 

96.1 

2.1 

98.8 

0.8 

98.5 

0.09 

99.1 


(6) Locate the candidate regions of the image awaiting detec¬ 
tion, obtain the feature vector using the foregoing methods, 
and input the computed feature vector into the trained SVM 
classifier. 

(7) Obtain the classification results and locate the pedestrian 
regions in the input image. 

6. Experimental results 

In this section, we report our experimental results with three 
different thermal infrared image databases from the OTCBVS 
Benchmark Dataset Collection [13]. The first and second test dat¬ 
abases are from the OSU Thermal Pedestrian Database [26] and 
the OSU Color-Thermal Database [27] separately. The first data¬ 
base contains 284 images with 155 pedestrians in total which 
were taken at an idle pedestrian intersection on the Ohio State 
University campus, while the second database contains 400 
images with 700 pedestrians which were taken at a busy pedes¬ 
trian intersection. The third test database is from the Terravic 
Motion IR Database provided by Roland Miezianko [13], taken 


in an outdoor environment, and containing 350 images with 
700 pedestrians. All the three databases are divided into train 
and test sets as usual. 

In order to evaluate the most efficient set of parameters Wi and 
w 2 listed in Eq. (1), we try different values for each parameter, both 
from 0.05 to 1 with the step 0.05, to get the flexible threshold T. We 
then use the different threshold T to segment the train infrared 
images. It is found that we can get the optimal segmentation re¬ 
sults when setting Wi and w 2 to 0.25 and 0.50 separately. Mean¬ 
while, we acquire the vertical and horizontal projection threshold 
T v and T h according to the statistics of the width/height ratios of 
376 pedestrians in the train sample set. In the experiments, the 
values of T v and T h are set to 0.66 and 0.30 separately. 

The experiments are divided into two groups. The first group 
compares the proposed wavelet entropy feature-based detection 
method with other well-known feature-based detection methods; 
and the second group compares the performance of different 
wavelet transforms. We use the True Pedestrians Rate (TPR) and 
the False Pedestrians Rate (FPR) to estimate the performance of 
the detection algorithm, where: 




a bed 

Fig. 11. Test images from database 3 using different wavelet transforms: (a) DWT, (b) DT CWT, (c) DD DWT, and (d) DD-DT CWT. 



Fig. 12. Test images from database 3 using different wavelet transforms (false alarm case): (a) DWT, (b) DT CWT, (c) DD DWT, and (d) DD-DT CWT. 
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„„„ true pedestrians count 

TPR =-- . - - . and 

pedestrians in total 

F p R _ false pedestrians count 

- non - pedestrians in total 

6.1. Comparisons between well-known features 

The well-known features used at present for pedestrian detec¬ 
tion in thermal infrared images are shape features (SHAPE) [10], 
histograms of oriented gradients (HOG) [11,12], and intensity dis¬ 
tribution-based inertia features (INERTIA) [9]. Three thermal infra¬ 
red image databases are tested with all these features and their 
potential hybrid features. Table 1 summarizes the experimental re¬ 
sults. Fig. 6 illustrates the performance comparisons of the differ¬ 
ent features based classifiers. From Table 1 and Fig. 6 we can see 
that the hybrid features using HOG + SHAPE + INERTIA is superior 
compared to other traditional features, while our proposed wavelet 
entropy (WE) feature-based classifier is superior compared to hy¬ 
brid features. The TPR of the WE feature is 1.02 times of the hybrid 
features, while the FPR of the WE feature is only one thirteenth of 
the hybrid features. Therefore our proposed feature has the best 
detection performance. Figs. 7 and 8 show some detection results 
for real thermal infrared images from test databases 2 and 3 sepa¬ 
rately (database 1 is ignored because of its simple background and 
the negligible differences between the experimental results). 

6.2. Comparisons between different wavelet transforms 

This section presents the comparisons between our proposed 
DD-DT CWT and other wavelet transforms, such as DWT, DT 
CWT and DD DWT. The wavelet entropy feature is used in all four 
wavelet transforms. Table 2 summarizes the experimental results. 
Fig. 8 shows the comparison of the performances of different wave¬ 
let transform based classifiers. We can also see that the proposed 
DD-DT CWT based classifier has the best detection performance 
from Table 2 and Figs. 9-12 show some detection results for real 
thermal infrared images from test databases 2 and 3 separately 
using the different wavelet transforms. 

7. Conclusions 

In this paper, a robust method for pedestrian detection in ther¬ 
mal infrared images has been proposed. The main characteristic of 
this method is the wavelet entropy feature-based classification 
algorithm using DD-DT CWT. The DD-DT CWT combines the DT 
CWT and DD DWT and has the advantages of shift invariance, 
directional selectivity, freedom from aliasing, and a near-continu¬ 
ous wavelet transform. The proposed wavelet entropy properly re¬ 
flects the energy distribution of the image in the frequency domain 
of the wavelet transform. The combination of DD-DT CWT and the 
wavelet entropy can thus describe the image features accurately. 
Experimental results have shown that our approach is very encour¬ 
aging. Future work will consider the temporal correlation and mo¬ 
tion cues for dealing with image sequences. Moreover, we will 
combine detection and tracking to improve the performance of 
our system. 
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